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Abstract 

Introduction: Several studies have reported the relationship between the 
STAT4 rs7574865G > T polymorphism as a susceptibility factor to ulcerative 
colitis (UC). However, the results have been controversial. Therefore, we 
conducted this meta-analysis to obtain the most reliable estimate of the 
association. 

Material and methods: PubMed, Embase and Web of Science databases 
were searched. Crude odds ratios (OR) with 95% confidence intervals (CI) 
were extracted and pooled to assess the strength of the association be- 
tween the STAT4 rs7574865G > T polymorphism and risk of UC. A total of 
five eligible studies including 1532 cases and 3786 controls based on the 
search criteria were involved in this meta-analysis. 

Results: We observed that the STAT4 rs7574865G > T polymorphism was 
significantly correlated with UC risk when all studies were pooled into the 
meta-analysis (the allele contrast model: OR = 1.13, 95% CI = 1.02-1.25; 
the heterozygote codominant model: OR = 1.22, 95% CI = 1.04-1.43; the 
dominant model: OR = 1.25, 95% CI = 1.07-1.45). In the stratified analysis 
by ethnicity, significant associations were observed in Spanish for the allele 
contrast model (OR = 1.20; 95% CI = 1.04-1.39), for the homozygote co- 
dominant model (OR = 1.57; 95% CI = 1.07-2.31), for the dominant model 
(OR = 1.20; 95% CI = 1.01-1.43), and for the recessive model (OR = 1.50; 
95% CI = 1.03-2.19). 

Conclusions: This meta-analysis suggests that the STAT4 rs7574865G > T 
polymorphism is a low-penetrant risk factor for UC, especially in Spanish. 
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Introduction 

Inflammatory bowel disease (IBD), including Crohn's disease (CD) 
and ulcerative colitis (UC), is a common health problem worldwide. 
It results from an interaction between genetic, environmental, inflec- 
tional, and other factors [1]. Preliminary findings suggest that genetic 
susceptibility plays a large part in the predisposition to UC [2, 3], and 
even the rare cases with Takayasu arteritis may display varietal genetic 
susceptibility [4, 5]. 

The signal transducer and activator of transcription-4 (STAT4) gene 
was found to be associated with multiple autoimmune diseases, such 
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as systemic lupus erythematosus, rheumatoid 
arthritis, type 1 diabetes, Wegener's granuloma- 
tosis, and ulcerative colitis [6]. Mutated allele T 
of the STAT4 rs7574865 SNP plays a key role in 
these autoimmune diseases [7, 8]. More recent- 
ly, several studies have assessed the relationship 
between the polymorphism of STAT4 G > T and 
the susceptibility to UC, but the results have been 
controversial [9-12]. Therefore, we conducted this 
meta-analysis to confirm these associations. 

Material and methods 

Identification and eligibility of relevant 
studies 

The search terms "STAR", "ulcerative colitis", 
"genotype", "polymorphism" and "variant" were 
employed to explore publications in PubMed, ISI 
Web of Knowledge and Embase databases for rel- 
evant reports (last search update October 2011). 
Only studies published in the English language 
were included. We did not define any minimum 
number of patients to be included for meta-anal- 
ysis. When multiple studies of the same patient 
population were identified, we included the pub- 
lished report with the largest sample size (Table I). 

Inclusion and exclusion criteria 

The following inclusion criteria were used to 
select literature for this analysis: (a) evaluation of 
the STAT4 rs7574865G > T polymorphism and UC, 
(b) only case-control studies were considered, (c) 
sufficient published data to estimate an odds ra- 
tio (OR) with 95% confidence interval (CI). Major 
exclusion criteria were: (1) no control population, 
(2) no available genotype frequency, and (3) dupli- 
cated studies. 

Data extraction 

For each study, the following data were collect- 
ed: first author's surname, year of publication, 
country of origin, ethnicity, criteria of enrolled pa- 
tients, genotyping method, total numbers of cases 
and controls as well as numbers of cases and con- 
trols with GG, GT and TT genotypes. The strength 
of the association between STAT4 rs7574865G > T 
polymorphism and UC risk was estimated using 
OR, with the corresponding 95% CI. Disagreement 
was resolved by discussion until a consensus was 
reached between the two authors. We did not de- 
fine any minimum number of patients for inclu- 
sion in our meta-analysis. 

Statistical analysis 

The risk of UC associated with STAT4 
rs7574865G > T was estimated for each study 
by OR with 95% CI. For five studies, we analysed 



the relationship for the allele contrast model (T 
vs. G). At the same time, due to lack the specific 
genotypes of Glas's literature reported by Glas, 
we estimated the association under other four 
different types of OR, namely the homozygote 
codominant model (TT vs. GG), the heterozygote 
codominant model (GT vs. GG), the dominant 
model (TT+GT vs. GG) and the recessive mod- 
el (TT vs. GT+GG). Hardy-Weinberg equilibrium 
(HWE) was tested by the x 2 test. The Q-statistic 
was used to investigate the degree of heteroge- 
neity between the trials, and a p-value of 0.10 
for the Q-test indicated a lack of heterogeneity 
among studies. We used the fixed-effects mod- 
el and the random-effects model based on the 
Mantel-Haenszel method [13] and the DerSimo- 
nian and Laird method [14], respectively, to com- 
bine values from each of the studies. A sensitivi- 
ty analysis was also performed by omitting each 
study in turn to identify potential outliers. All of 
the statistical analyses were performed with Re- 
view Manage version 4.3 and STATA version 12.0 
using two-sided p-values. 

Results 

Study characteristics 

We obtained 12 studies about the association 
between STAT4 rs7574865 polymorphism. Fol- 
lowing the above inclusion and exclusion criteria, 
5 publications were included in the final meta- 
analysis [9-12]. Characteristics of studies focusing 
on STAT4 rs7574865G > Tare summarized in Table II. 
Genotypes and separate p values for STAT4 
rs7574865 polymorphism are listed in Table III. 

Meta-analysis results 

Table III lists the main results of pooled ORs 
for STAT4 rs7574865G > T polymorphism and UC 
risk. We observed that the STAT4 rs7574865G > T 
polymorphism was significantly correlated with 
UC risk when all studies were pooled into the me- 
ta-analysis (the allele contrast model: OR = 1.13, 
95% CI = 1.02-1.25, p = 0.224 for heterogenei- 
ty, Figure 1; the heterozygote codominant model: 
OR = 1.22, 95% CI = 1.04-1.43, p = 0.49 for het- 
erogeneity, Figure 2; the dominant model: OR = 
1.25, 95% CI = 1.07-1.45, p = 0.54 for heterogene- 
ity, Figure 3). In the analysis stratified by ethnicity, 
significant associations were observed in Spanish 
for the allele contrast model (OR = 1.20; 95% CI 
= 1.04-1.39; p = 0.42 for heterogeneity), for 
the homozygote codominant model (OR = 1.57; 
95% CI = 1.07-2.31; p = 0.69 for heterogeneity), 
for the dominant model (OR = 1.20; 95% CI = 
1.01-1.43; p = 0.40 for heterogeneity), and for the 
recessive model (OR = 1.50; 95% CI = 1.03-2.19; 
p = 0.78 for heterogeneity). 
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Table I. Characteristics of UC studies included in the meta-analysis 



First author/published year 


Study population 
of cases 


Method 


Source 


Cases 


Controls 


HWE 
of controls 


Bouzid [9] 2009 




North Africa 


SNP genotyping 
assays 


HB 


68 


162 




0.494 


Diaz-Gallo [8] 2010 




Spanish 


SNP genotyping 
assays 


HB 


402 


1296 




0.888 


Moon [7] 2010 




Korean 


SNP genotyping 
assays 


HB 


246 


229 




0.934 


Martinez [10] 2008 




Spanish 


SNP genotyping 
assays 


HB 


352 


716 




0.968 


Glas [11] 2010 




Caucasian 


SNP genotyping 
assays 


HB 


464 


1383 




0.215 


Table II. Genotypes and p-values of STAT4 rs7574865G > T polymorphism included in the meta analysis 


First author 




Cases 








Controls 








GG 


GT TT 


G T 


GG 


GT 


TT 


G 


T 


Bouzid 


39 


25 4 


103 33 


113 


46 


3 


272 


52 


Diaz-Gallo 


241 


137 24 


619 185 


813 


428 


55 


2054 


538 


Moon 


97 


131 18 


325 167 


106 


100 


23 


312 


146 


Martinez 


209 


123 20 


541 163 


470 


220 


26 


1160 


272 


Glas 






729 199 








2171 


595 



Table III. Results of meta-analysis for STAT4 rs7574865 polymorphisms and UC risk 



Study groups 


Allele contrast model 


Homozygote codominant 






OR (95% CI) 


Value of p 


OR (95% CI) 


Value of p 




Total 


1.13 (1.02-1.25) 


0.006 


1.39 (1.00-1.92) 


0.05 




Ethnicity 


Spanish 


1.20 (1.04-1.39) 


0.01 


1.57 (1.07-2.31) 


0.02 




Korean 


1.10 (0.84-1.44) 




0.86 (0.44-1.68) 






North Africa 


1.33 (0.82-2.15) 




3.86 (0.83-18.03) 






Study groups 


Heterozygote codominant 


Dominant model 


Recessive mode 






OR (95% CI) 


Value 
of p 


OR (95% CI) Value 
of p 


OR (95% CI) 


Value 
of p 


Total 


1.22 (1.04-1.43) 


0.01 


1.25 (1.07-1.45) 0.004 


1.30 (0.82-2.08) 


0.26 


Ethnicity 


Spanish 


1.15 (0.96-1.38) 


0.12 


1.20 (1.01-1.43) 0.04 


1.50 (1.03-2.19) 


0.04 


Korean 


1.43 (0.98-2.09) 




1.32 (0.92-1.91) 


0.71 (0.37-1.35) 




North Africa 


1.57 (0.86-2.89) 




1.71 (0.95-3.08) 


3.31 (0.72-15.22) 





Publication bias 

We performed Begg's funnel plot and Egger's 
test to assess the publication bias of the litera- 
ture. The results did not show any evidence of 
publication bias in all the comparisons. We pres- 
ent funnel plot for ORs of C vs. T in Figure 4. Also, 
the results of Egger's test still did not suggest any 



evidence of publication bias (p = 0.679 for T vs. G; 
p = 0.134 for GT vs. GG; p = 0.104 for TC vs. TT). 

Discussion 

The Stat signalling pathways play important 
roles in carcinogenesis and immunopathology of 
many malignancies and immunologic diseases. 
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Study ID 



OR (95% CI) % Weight 



Bouzid (2009) 
Diaz-Gallo (2010) 
Moon (2010) 
Martinez (2008) 
Glas (2010) 

Overall (/ 2 = 29.7%, p = 0.224) 



1.68 (1.03,2.74) 3.37 



1.14 (0.94,1.38) 28.35 



1.10 (0.84, 1.44) 14.44 



1.28 (1.03, 1.60) 19.91 



1.00 (0.83, 1.19) 33.93 



1.13 (1.02, 1.25) 100.00 



0.5 



Figure 1. Forest plot for the association between the STAT4 rs7574865G > T polymorphism and ulcerative colitis 
risk in overall studies (for G vs. T allele contrast model). The study is shown by the point estimate of the OR (the 
size of the square is proportional to the weight of each study) and 95% CI for OR; pooled OR and 95% CI have been 
appropriately derived from the fixed effects model 



Study ID 



OR (95% CI) % Weight 



Martinez (2008) 



Bouzid (2009) 



Diaz-Gallo (2010) 



Moon (2010) 



Overall (/ 2 = 0.0%, p = 0.494) 



1.26 (0.96, 1.65) 32.34 



1.57 (0.86, 2.89) 



1.08 (0.85, 1.37) 



5.78 



45.80 




1.43 (0.98, 2.09) 16.07 



1.22 (1.04, 1.43) 100.00 



0.9 



Figure 2. Forest plot for the association between the STAT4 rs7574865G > T polymorphism and ulcerative colitis 
risk in overall studies (for the heterozygote codominant model, fixed effects model) 



Recent studies have implicated their members' 
possible involvement in the pathogenesis of IBD. 
For example, the expression and/or activation of 
interleukin 12 (IL-12) and Stat4, including phos- 
phorylated Stat4, suggest that proinflammato- 
ry IL-12/Stat4 signalling or IL-23/Stat_4 are likely 
candidate pathways involved in the inflamma- 
tory pathology in UC and colorectal cancer (CRC) 
[15-17]. Moreover, it is generally recognized that 
there is an increased risk of CRC in patients with 
UC, and the overall prevalence of manifest CRC 
in patients with UC is unacceptably high. So, this 
understanding may have an important impact on 
patients with certain risk factors, such as STAT4, 
etc [18]. Some researchers have examined the as- 
sociation of STAT4 rs7574865G > T polymorphism 



with UC risk and treatments. Based on previous 
studies reporting increased sensitivity to IFN-a in 
lupus patients carrying the risk variant of STAT4 
might contribute to increased mucosal inflamma- 
tion in IBD patients and to the response to immu- 
nosuppressive and immunomodulatory therapies 
[19], and a significant relationship was observed 
in several but not all studies [20, 21]. 

The aim of our study was to demonstrate the 
role of STAT4 rs7574865G > T polymorphism in the 
relationship with UC risk using a meta-analysis. 
Interestingly, we observed that STAT4 rs7574865G 
> T was significantly correlated with UC risk and 
the T allele of the STAT4 rs7574865G > T variant 
was a low-penetrant risk factor for UC. Simultane- 
ously there was evidence to indicate that STAT4 
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Study ID 



OR (95% CI) % Weight 



Martinez (2008) 



Bouzid (2009) 



Diaz-Gallo (2010) 



Moon (2010) 



Overall (/ 2 = 0.0%, p = 0.543) 




1.31 (1.01, 1.70) 32.07 



1.71 (0.95,3.08) 5.53 



1.12 (0.89, 1.41) 45.67 



1.32 (0.92, 1.91) 16.73 



1.25 (1.07, 1.45) 100.00 



I 

0.9 1 3 
Figure 3. Forest plot for the association between the STAT4 rs7574865G > T polymorphism and ulcerative colitis 
risk in overall studies (for the dominant model, fixed effects model) 



rs7574865G > T polymorphism was associated with 
increased risk of UC in Spanish and mixed ethni- 
cities. 

Of course, we must note that our study pop- 
ulation was small, so it is still necessary to con- 
duct larger sample studies considering gene-gene 
and gene-environment interactions, and using 
standardized unbiased genotyping methods, ho- 
mogeneous UC patients, and sufficiently matched 
controls. In addition, we also should consider the 
relationship between treatment response and 
STAT4 polymorphisms in various types of UC in 
future work for a comprehensive evaluation of 
the relevance of this polymorphism for UC clinical 
characteristics, for example, whether it is similar 
to TLR2 polymorphisms that can predict the inci- 
dence of steroid-dependent UC and the response 
to azathioprine treatment, etc [12, 22, 23]. 

Recently, genome-wide association studies 
(GWAS) have helped scientists understand the in- 
heritance patterns of disorders on a global scale. 
In particular, this may lead to an enormous boost 
to the identification of susceptibility genes for 
several diseases, including inflammatory bowel 
disease, but there are some insurmountable prob- 
lems, the most critical being that some high-risk 
sites for low frequency are masked by low-risk 
sites for high frequency, so this may result in the 
analysis not being comprehensive, with the loss of 
a lot of meaningful SNP sites. For example, Rivas 
et al. have identified more than 1,000 suscepti- 
bility loci for IBD through GWAS; further analysis 
found that only four additional independent risk 
factors, NOD2, two additional protective variants 
in IL23R, and a protective splice variant in CARD9, 
are associated with IBD [21], while some of the 
classic sites, such as HLA, cytokine and STAT4 gene 
polymorphisms [24, 25], cannot be identified. 



O 0.05 



o 0.10- 



-a 0.15 



ft 0.20- 



0.25 




-0.4 



-0.2 



0 0.2 
Log OR 



0.4 



0.6 



Figure 4. Funnel plot analysis for odds ratios of G 
allele compared with T allele overall 

In conclusions, our study shows a genetic as- 
sociation between STAT4 rs7574865G > T vari- 
ant and increased susceptibility to UC. However, 
it needs to be further evaluated in larger sample 
collections in diverse ethnic populations and more 
types of UC along with clinical characteristics, to 
confirm our findings. 
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